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CONSTRAINTS ON THE COSMIC RAY DIFFUSION COEFFICIENT 
IN THE W28 REGION FROM GAMMA-RAY OBSERVATIONS 

S. Gabici 1 , S. Casanova 2 ' 3 , F. A. Aharonian 2 ' 4 and G. Rowell 5 

Abstract. GeV and TeV gamma rays have been detected from the supernova remnant W28 and its 
surroundings. Such emission correlates quite well with the position of dense and massive molecular clouds 
and thus it is often interpreted as the result of hadronic cosmic ray interactions in the dense gas. Constraints 
on the cosmic ray diffusion coefficient in the region can be obtained, under the assumption that the cosmic 
rays responsible for the gamma ray emission have been accelerated in the past at the supernova remnant 
shock, and subsequently escaped in the surrounding medium. In this scenario, gamma ray observations can 
be explained only if the diffusion coefficient in the region surrounding the supernova remnant is significantly 
suppressed with respect to the average galactic one. 
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Cosmic ray (CR) spallation measurements allow us to infer the average residence time t res of CRs in the 
Galaxy. Then, if h is the distance that a CR has to move away from its source before escaping the Galaxy 
(i.e. the Galaxy's thickness), it is possible to estimate the CR diffusion coefficient asflsi h 2 /t res . A thorough 
comparison between spallation data and CR transport models gives: D ga \ w 10 28 (£'/10 GeV) 5 cm 2 s _1 , where 
E is the CR particle energy and S ~ 0.3 -j- 0.7 (Berezinskii et al. 1990[ ). However, this has to be intended as 
the average CR diffusion coefficient in the Galaxy, and local variations (both in time and space) of the diffusion 
coefficient might exist. 

In particular, the diffusion coefficient might be suppressed close to CR sources. This is because CRs can excite 
magnetic turbulence while streaming away from their acceleration site. This would enhance the scattering rate 
of CR themselves and thus reduce the diffusion coefficient ( Kulsrud fc Pierce|1969 1. The problem of estimating, 
on theoretical grounds, the diffusion coefficient around CR sources is far from being solved, mainly because of 
its intrinsic non-linearity and because various mechanisms might damp the CR generated waves and thus affect 
the way in which CRs diffuse (e.g. |Farmer k Goldreich|[2004| |Ptuskin et aL|2008 ). 

Gamma ray observations can provide us with constraints on the diffusion coefficient close to CR sources. 
After escaping from their sources, CRs undergo hadronic interactions with the surrounding gas and produce 
gamma rays. The characteristics of such radiation (i.e. its spectrum and intensity as a function of the time 
elapsed since CRs escaped the source) depend on the value of the diffusion coefficient that can thus be con- 
strained, if a reliable modeling for CR acceleration at the source is available. The presence of massive Molecular 
Clouds (MCs) close to the source would enhance the gamma ray emission, making its detection more probable. 

The study of gamma ray emission from runaway CRs has been studied by Aharonian & Atoyan ( 1996 ) for 



the case of a generic CR source, while the specific (and most popular) situation of CR acceleration at SuperNova 
Remnants (SNRs) has been described in a series of recent papers ( Gabici fc Aharonian|2007 Gabici et al.|2 009; 
Casanova et al. 20101. The study of such radiation is of great importance not only in order to reach a better 



understanding of the way in which CRs diffuse in the interstellar magnetic field, but also because its detection 
can provide an indirect way to identify the sources of galactic CRs. 
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W28 is an old SNR in its radiative phase of evolution, located in a region rich of dense molecular gas with 
average density > 5 cm -3 . At an estimated distance of ~ 2 kpc the SNR shock radius is ~ 12 pc and its 
velocity ~ 80 km/s (e.g. Rho & Borkowski 2002). By using the dynamical model by|Ciom et al. (1988) and 



assuming that the mass of the supernova ejecta is ~ 1.4 M , it is possible to infer the supernova explosion 
energy (E SN ~ 0.4 x 10 51 erg), initial velocity (~ 5500 km/s), and age (t age ~ 4.4 x 10 4 yr). 

Gamma ray emission has been detected from the region surrounding W28 both at TeV (Aharonian et al. 
2008| and GeV energies ( |Abdo et al.poTo) |Giuliani et al.||2010| , by HESS, FERMI, and AGILE, respectively. 
The TeV emission correlates quite well with the position of three massive molecular clouds, one of which is 
interacting with the north-eastern part of the shell (and corresponds to the TeV source HESS J1801-233), and 
the other two being located to the south of the SNR (TeV sources HESS J1800-240 A and B) . The masses 
of these clouds can be estimated from CO measurements and result in s» 5, 6, and 4 x 10 4 M Q , respectively, 
and their projected distances from the centre of the SNR are ps 12, 20, and 20 pc, respectively (Aharonian et 



al. 20081. The GeV emission roughly mimics the TeV one, except for the fact that no significant emission is 



detected at the position of HESS J1800-240 A. 

In this paper, we investigate the possibility that the gamma ray emission from the W28 region could be the 
result of hadronic interactions of CRs that have been accelerated in the past at the SNR shock and then escaped 
in the surrounding medium Q To do so, we follow the approach described in Gabici et al. (2009) (hereafter 
GAC2009), which we briefly summarize below. 

For each particle energy E we solve the diffusion equation for CRs escaping the SNR. For simplicity we 
treat the SNR as a point like source of CRs and we consider an isotropic and homogeneous diffusion coefficient 
D = xDgai oc E - 5 (see Eq. 5 in GAC2009). Here, \ represents possible deviations with respect to the average 
CR diffusion coefficient in the Galaxy. The solution of the diffusion equation gives the spatial distribution of CRs 
around the source fcR, which is roughly constant up to a distance equal to the diffusion radius Rd — y/i D tdiff, 
and given by fcR oc tjEsn / 'R% where r\ is the fraction of the supernova explosion energy converted into CRs, 
and tdiff is the time elapsed since CRs with energy E escaped the SNR (see eq. 6 in GAC2009). For distances 
much larger than Rd the CR spatial distribution falls like fcR oc exp(—(R/Rd) 2 ) (see eq. 6 in GAC2009). 
Following GAC2009, we assume that CRs with energy 5 PeV (1 GeV) escape the SNR at the beginning (end) 
of the Sedov phase, at a time ~ 250 yr (~ 1.2 x 10 4 yr) after the explosion, and that the time integrated CR 
spectrum injected in the interstellar medium is oc E~ 2 . In this scenario, particles with lower and lower energies 
are released gradually in the interstellar medium ( Ptuskin fc Zirakashvili|2005 1 , and we parametrize the escape 
time as: t esc oc E~ a which, during the Sedov phase, can also be written as R s oc E~ 2a / 5 , where R s is the shock 
radius at time t esc and a ~ 4. From this it follows that the assumption of point like CR source is good for high 
energy CRs only (~ TeV or above), when R s is small, but it becomes a rough approximation at significantly 
lower energies. This is because low energy particles are believed to be released later in time, when the SNR 
shock radius is large (i.e. non negligible when compared to Rd). 

We now provide a simplified argument to show how we can attempt to constrain the diffusion coefficient by 
using the TeV gamma ray observations of the MCs in the W28 region. The time elapsed since CRs with a given 
energy escaped the SNR can be written as: tdiff = t age — t esc . However, for CRs with energies above 1 TeV 
(the ones responsible for the emission detected by HESS) we may assume t esc « t age (i.e. high energy CRs 
are released when the SNR is much younger than it is now) and thus tdiff ~ t age . Thus, the diffusion radius 
reduces to Rd ~ ^/4 D t age . We recall that within the diffusion radius the spatial distribution of CRs, fcR, is 
roughly constant, and proportional to rjEsN/Rd- O n the other hand, the observed gamma ray flux from each 
one of the MCs is: F 1 oc fcR^ci/d 2 , where M c \ is the mass of the MC and d is the distance of the system. 
Note that in this expression F 1 is calculated at a photon energy E 1: while fcR is calculated at a CR energy 
Ecr ~ 10 x E~j, to account for the inelasticity of proton-proton interactions. By using the definitions of fcR 
and Rd we can finally write the approximate equation, valid within a distance Rd from the SNR: 



V E S n ( Md 
(x D gal t age )^ \ d? 



Estimates can be obtained for all the physical quantities in the equation except for the CR acceleration 
efficiency r\ and the local diffusion coefficient xDgai ■ By fitting the TeV data we can thus attempt to constrain, 
within the uncertainties given by the errors on the other measured quantities (namely, Esn, t age , M c i, and d) 



"This scenario has been described in a number of recent papers | |Fujita et al.|2009| [Ohira et al.|2010||Li fc Chen|2010| 



Cosmic ray diffusion in the W28 region 
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Fig. 1. Simultaneous fit to the three TeV sources detected by HESS in the W28 region. Gamma ray spectra have been 
calculated by using the parameterizations by Kamae et al. (20061, where a multiplicative factor of 1.5 has been applied 
to account for the contribution to the emission from nuclei heavier than H both in CRs and in the interstellar medium. 



and by the assumptions made (e.g. the CR injection spectrum is assumed to be E~ 2 ), a combination of these 
two parameters (namely rj/x 3 ^ 2 )- The fact that the MCs have to be located within a distance Rd from the SNR 
can be verified a posteriori. Given all the uncertainties above, our results have to be interpreted as a proof of 
concept of the fact that gamma ray observations of SNR/MC associations can serve as tools to estimate the 
CR diffusion coefficient. More detection of SNR/MC associations are needed in order to check whether the 
scenario described in this paper applies to a whole class of objects and not only to a test-case as W28. Future 
observations from the Cherenkov Telescope Array will most likely solve this issue. 

Fig. [I] shows a fit to the HESS data for the three massive MCs in the W28 region. A simultaneous fit to 
all the three MCs is obtained by setting Tj/x 3 ^ 2 ~ 20, which implies that the diffusion coefficient (normalized 
to the average galactic one) x nas to be much smaller than 1 for any reasonable value of rj < 1. For example, 
an acceleration efficiency rj = 30% corresponds to a CR diffusion coefficient of x — 0.06, which in turn gives a 
diffusion distance for TeV particles of Rd ~ 60 pc. This means that the results in Fig. [T] are valid if the physical 
(not projected) distances between the MCs and the SNRs do not significantly exceed Rd- Small values of the 
diffusion coefficient have been also proposed by Giuliani et al. (2010); Fujita et al.| ( |2009~ ); Li & Chen (2010). 
Note that, since we are considering gamma rays in a quite narrow (about one order of magnitude) energy band 
around ~ 1 TeV, we are actually constraining the diffusion coefficient of CRs with energy ss 10 TeV, and we 
cannot say much about the energy dependence of the diffusion coefficient on a broad energy interval. 

In principle, observations by FERMI and AGILE might be used to constrain the diffusion coefficient down to 
GeV particle energies. However, in this energy band the uncertainties are more severe because of the following 
reasons: i) low energy CRs are believed to be released late in time, when the SNR shock is large, and thus the 
assumption of point-like source is probably not well justified (see Ohira et al.||2010 for a model that takes into 
account the finite size of the SNR) ; ii) for the same reason, we can no longer assume that tdiff ~ t age , as we did 
for high energy CRs. In other words, we need to know the exact time at which CRs with a given energy escape 
the SNR. Though some promising theoretical studies exist (e.g. Ptuskin fc Zirakashvili|20'05 1 , our knowledge of 
the escape time of CRs from SNRs is still quite limited. 

Fig. [2] shows a fit to the broad band gamma ray spectrum measured from FERMI and HESS. The three 
panels refers to (left to right) the sources HESS J1801-233, HESS J1800-240 A and B, respectively. Dashed lines 
represent the contribution to the emission from CRs that escaped from W28, dotted lines the contribution from 
background CRs, and solid lines the total emission. Since FERMI data refers to the emission after background 
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Fig. 2. Broad band fit to the gamma ray emission detected by FERMI and HESS from the sources HESS J1801-233, 
HESS J1800-240 A and B (left to right). Dashed lines represent the contribution to the gamma ray emission from 
CRs that escaped W28, dotted lines show the contribution from the CR galactic background, and solid lines the total 
emission. Distances to the SNR centre are 12, 65, and 32 pc (left to right). FERMI and HESS data points are plotted 
in black. No GeV emission has been detected from HESS J1800-240 A. 

subtraction, dashed lines have to be compared with data points. The (often non-trivial) background subtraction 
issue might add another source of uncertainty in the comparison between data and predictions. An acceleration 
efficiency r\ = 30% and a diffusion coefficient \ — 0-06 have been assumed, while the distance from the SNR 
centre is assumed to be (left to right) 12, 65, and 32 pc. Keeping in mind all the above mentioned caveats, it 
is encouraging to see that a qualitative agreement exists between data and predictions also in the GeV band. 

We investigated the possibility that the gamma ray emission detected from the MCs in the region of the 
SNR W28 is produced by CRs that escaped the SNR. This interpretation requires the CR diffusion coefficient 
in that region to be significantly suppressed with respect to the average galactic one. Such suppression might 
be the result of an enhancement in the magnetic turbulence due to CR streaming away from the SNR. 

SG thanks M. Walmsley for many interesting discussions. The research leading to these results has received funding from the 
European Union [FP7/2007-2013] under grant agreement n° 256464 and 221093. 
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